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ABSTRACT 

Aims. We investigate the impact of polarised foreground emission on the performances of future CMB experiments aiming the detec- 
tion of primordial tensor fluctuations in the early universe. In particular, we study the accuracy that can be achieved in measuring the 
tensor-to-scalar ratio r in presence of foregrounds. 

Methods. We design a component separation pipeline, based on the Smica method, aimed at estimating r and the foreground con- 
tamination from the data with no prior assumption on the frequency dependence or spatial distribution of the foregrounds. We derive 
error bars accounting for the uncertainty on foreground contribution. We use the current knowledge of galactic and extra-galactic fore- 
grounds as implemented in the Planck Sky Model (PSM), to build simulations of the sky emission. We apply the method to simulated 
observations of this modelled sky emission, for various experimental setups. 

Results. Our method, with Planck data, permits us to detect r = 0. 1 from B-modes only at more than 3cr. With a future dedicated 
space experiment, as EPIC, we can measure r = 0.001 at ~ 6<r for the most ambitious mission designs. Most of the sensitivity to 
r comes from scales 20 < £ < 150 for high r values, shifting to lower f s for progressively smaller r. This shows that large scale 
foreground emission doesn't prevent a proper measurement of the reionisation bump for full sky experiment. We also investigate the 
observation of a small but clean part of the sky. We show that diffuse foregrounds remain a concern for a sensitive ground-based 
experiment with a limited frequency coverage when measuring r < 0. 1 . Using the Planck data as additional frequency channels to 
constrain the foregrounds in such ground-based observations reduces the error by a factor two but does not allow to detect r = 0.01. 
An alternate strategy, based on a deep field space mission with a wide frequency coverage, would allow us to deal with diffuse fore- 
grounds efficiently, but is in return quite sensitive to lensing contamination. In the contrary, we show that all-sky missions are nearly 
insensitive to small scale contamination (point sources and lensing) if the statistical contribution of such foregrounds can be modelled 
accurately. Our results do not significantly depend on the overall level and frequency dependence of the diffused foreground model, 
when varied within the limits allowed by current observations. 

Key words, cosmic microwave background - cosmological parameters - Cosmology: observations 
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<T^> 1. Introduction however, the only one. Foreground emissions from the galac- 
tic interstellar medium (ISM) and from extra-galactic objects 

# £ 'After the success of the WMAP space mission in map- (galaxies and clusters of galaxies) superimpose to the CMB. 

k> i ping the Cosmic Microwave Background (CMB) temperature Most foregrounds are expected to emit polarised light, with a 

•jjj • anisotropies, much attention now turns towards the challenge polarisation fraction typically comparable, or larger, than that 

C3 ; of measuring CMB polarisation, in particular pseudo-scalar po- f the CMB. Component separation (disentangling CMB emis- 

larisation modes (the B-modes) of primordial origin. These s i on f rom a n these foregrounds) is needed to extract cosmologi- 

B-modes offer one of the best options to constrain infla- ca i information from observed frequency maps. The situation is 

tionary models flSeljak & Z aldarriagalfl997t iHu & White 1997; particularly severe for the B-modes of CMB polarisation, which 



Kamionkowski et all 119971: |Kamionkowski & Kosowsk^|1998|; w m be, if measurable, sub-dominant at every scale and every 
Baumann & Peirisll2008l) . frequency. 



First polarisation measurements have already be e n ob- T h e main objective of this paper is to evaluate the accuracy 
tained by a number of instruments ( Kovac et afl |2002t w ith which various upcoming or planned experiments can mea- 
ISievers & CBI Collaboration||2005|;|Page et al.||2007|), but no de- sure r in prese nce of foregrounds . This problem has been ad- 
tection of B-modes has been claimed yet. While several ground- dressed before: iTucci et alj d|005) investigate the lower bound 
based and balloon-borne experiments are already operational, for r that can be achieved considering a simple foreground clean- 
or in construction, no CMB-dedicated space-mission is planned ing technique, based on the extrapolation of foreground tern- 
after Planck at the present time: whether there should be one for p lates and subtraction from a channel dedicated to CMB mea- 
CMB B-modes, and how it should be designed, are still open sur ement: IVerde et alJ (120061) assume foreground residuals at a 
questions. known level in a cleaned map, treat them as additional Gaussian 
As CMB polarisation anisotropies are expected to be signifi- noise, and compute t he error on r due to such excess noise; 
cantly smaller than temperature anisotropies (a few per cent at lAmblard et al.l d2007l) investigate how best to select the fre- 
most), improving detector sensitivities is the first major chal- quency bands of an instrument, and how to distribute a fixed 
lenge towards measuring CMB polarisation B-modes. It is not, number of detectors among them, to maximally reject galactic 
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Fig. 1. Respective emission levels of the various components as predicted by the PSM. Left: predicted power spectra of the various 
components at 100 GHz, compar ed to CMB and lensi ng level for standard cosmology and various values of r (r = 0.07, and other 
cosmological parameters follow Dunklev et al. (2008a)). The power spectra of diffuse galactic foregrounds are computed using 
the cleanest 55% of the polarised sky. The power spectrum from residual point sources is computed assuming that all sources 
brighter than 500 mJy (in temperature) in one of the Planck channels have been cut out. Right: typical frequency-dependence of the 
contributions to B-type polarisation of CMB, synchrotron and dust, at 1 degree resolution. The dashed lines correspond to the mean 
level of fluctuation as computed outside the mask used for the power spectra shown in the right panel. 



foreground contamination. This latter analysis is based on an 
Intern al Linear Combinatio n cleaning technique similar to the 
one of Tegmar k et al.l d2003l) on WMAP temperature anisotropy 
data. The two last studies assume somehow that the residual con- 
tamination level is perfectly known - an information which is 
used to derive error bars on r. 

In this paper, we relax this assumption and propose a method 
to estimate the uncertainty on residual contamination from the 
data themselves, as would be the case for real data analysis. We 
test our method on semi-realistic simulated data sets, including 
CMB and realistic foreground emission, as well as simple instru- 
mental noise. We study a variety of experimental designs and 
foreground mixtures. 

This paper is organised as follows: the next section (Sect. |2]i 
deals with polarised foregrounds and presents the galactic emis- 
sion model used in this work. In section[3] we propose a method, 
using the most recent version of the Smica component separa- 
tion framework (ICardoso et al.ll2008l) . to provide measurements 
of the tensor to scalar ratio in presence of foregrounds. In sec- 
tion!?] we present the results obtained by applying the method to 
various experimental designs. Section[5]discusses the reliability 
of the method (and of our conclusions) against various issues, in 
particular modelling uncertainty. Main results are summarised in 
section|6] 

2. Modelling polarised sky emission 

Several processes contribute to the total sky emission in the fre- 
quency range of interest for CMB observation (typically be- 
tween 30 and 300 GHz). Foreground emission arises from the 
galactic interstellar medium (ISM), from extra-galactic objects, 
and from distortions of the CMB itself through its interaction 
with structures in the nearby universe. Although the physical 
processes involved and the general emission mechanisms are 
mostly understood, specifics of these polarised emissions in the 



millimetre range remain poorly known as few actual observa- 
tions, on a significant enough part of the sky, have been made. 

Diffuse emission from the ISM arises through synchrotron 
emission from energetic electrons, through free-free emission, 
and through grey-body emission of a population of dust grains. 
Small spinning dust grains with a dipole el ectric moment may 
also e mit significantly in the radio domain dDraine & Lazarianl 
1998). Among those processes, dust and synchrotron emissions 
are thought to be significantly polarised. Galactic emission also 
includes contributions from compact regions such as supernovae 
remnants and molecular clouds, which have specific emission 
properties. 

Extra-galactic objects emit via a number of different mecha- 
nisms, each of them having its own spectral energy distribution 
and polarisation properties. 

Finally, the CMB polarisation spectra are modified by the in- 
teractions of the CMB photons on their way from the last scatter- 
ing surface. Reionisation, in particular, re-injects power in polar- 
isation on large scales by late-time scattering of CMB photons. 
This produces a distinctive feature, the reionisation bump, in the 
CMB B-mode spectrum at low t. Other interactions with the lat- 
ter universe, and in particular lensing, contribute to hinder the 
measurement of the primordial signal. The lensing effect is par- 
ticularly important on smaller scales as it converts a part of the 
dominant E-mode power into B-mode. 

In the following, we review the identified polarisation pro- 
cesses and detail the model used for the present work, with a 
special emphasis on B-modes. We also discuss main sources of 
uncertainty in the model, as a basis for evaluating their impact 
on the conclusions of this paper. 

Our simulations are based on the Planck Sky Model (PSM), 
a sky emission simulation tool developed by the Planck col- 
laboration for pre-launch preparation of Planck data analysis 
dDelabrouille et all 12009). Figure Q] gives an overview of fore- 
grounds as included in our baseline model. Diffuse galactic 
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emission from synchrotron and dust dominates at all frequen- 
cies and all scales, with a minimum (relative to CMB) between 
60 and 80 GHz, depending on the galactic cut. Contaminations 
by lensing and a point source background are lower than primor- 
dial CMB for r > 0.01 and for I < 100, but should clearly be 
taken into account in attempts to measure r < 0.01. 

2.1. Synchrotron 

Cosmic ray electrons spiralling in the galactic magnetic 
field produce highly pola rised synchrotron emission (e.g. 
Rybicki & Lightmar]dl979j)). This is the dominant contaminant 
of the polarised CMB signal at low frequency (< 80 GHz), as 
can be seen in the right panel of Fig. [TJ In the frequency range 
of interest for CMB observations, measurements of this emis- 
sion have b een provided, both in temperature and polarisation, 
by WMAP dPage et alJl2007t iGold et al.ll2008l) . The intensity of 
the synchrotron emission depends on the cosmic ray density n e , 
and on the strength of the magnetic field perpendicularly to the 
line of sight. Its frequency scaling and its intrinsic polarisation 
fraction f s depend on the energy distribution of the cosmic rays. 

2.1.1. Synchrotron emission law 

For electron density following a power law of index p, n e (E) oc 
E~ p , the synchrotron frequency dependence is also a power law, 
of index /3 S = -(p + 3)/2: 



S(v) = S(v )(v/v o y* 



(1) 



where the spectral index, is equal to -3 for a typical value 
p = 3. 

The synchrotron spectral index depends significantly on cos- 
mic ray properties. It varies with the direction of the sky, and 
possib ly, with the frequency of observation (see e.g. lStrong et al.l 
(2007|) for a review of propagation and interaction processes of 
cosmic rays in the galaxy). 

For a multi-channel experiment, the consequence of this is 
a decrease of the coherence of the synchrotron emission across 
channels, i.e. the correlation between the synchrotron emission 
in the various frequency bands of observation will be below 
unity. 

Observational constraints have been put on the synchrotron 
emission law. A template of s ynchrotron emission intensity at 
408 MHz has been provided bv lkaslam et~ai](ll982l). Combining 
this map with sky surveys at 1.4 GHz dReich & Reich|[ l986) 
and 2.3 GHz d Jonas etal.1 1 1998b . iGiardino et ail d2002b and 
iPlatania et alj d2003 ) have derived nearly full sk y spectral in- 
dex m aps. Using the measurement from WMAP, iBennett et alj 
( 2003) derived the spectral index between 408 MHz and 23 GHz. 
Compared to the former results, it showed a significant steep- 
ening toward /3 S = -3 around 20 GHz, and a strong galac- 
tic plane feature with flatter spectral index. This feature was 
first interpreted as a flatter cosmic ray distribution in star form- 
ing regions. Recently, however, taking into account the pres- 
ence, at 23 GHz, of additional contribution from a possi- 
ble anomalous emission correlate d with the dust column den- 
sity, iMiville-Deschenes et al.l d2008b found no such pronounced 
galactic feature, in better agreement with lower frequency re- 
sults. The spectral index map obtained in this way is consistent 
with jS s = -3 ± 0.06. 

There is, hence, still significant uncertainty on the exact vari- 
ability of the synchrotron spectral index, and in the amplitude of 
the steepening if any. 



2.1.2. Synchrotron polarisation 

If the electron density follows a power law of index p, the syn- 
chrotron polarisation fraction reads: 



/, = 3(p + l)/(3p + 7) 



(2) 



For p = 3, we get f s - 0.75, a polarisation fraction which 
varies slowly for small variations of p. Consequently, the intrin- 
sic synchrotron polarisation fraction should be close to constant 
on the sky. However, geometric depolarisation arises due to vari- 
ations of the polarisation angle along the line of sight, partial 
cancellation of polarisation occurring for superposition of emis- 
sion with orthogonal polarisation directions. Current measure- 
ments show variations of the observed polarisation value from 
about 10% near the gala ctic plane, to 30-50 % at intermediate to 
high galactic latitudes dMacellari et al.ll2008h . 

2.1 .3. Our model of synchrotron 

In summary, the B-mode intensity of the synchrotron emission is 
modulated by the density of cosmic rays, the slope of their spec- 
tra, the intensity of the magnetic field, its orientation, and the 
coherence of the orientation along the line of sight. This makes 
the amplitude and frequency scaling of the polarised synchrotron 
signal dependant on the sky position in a rather complex way. 

For the purpose of the pr esent work, we mostly follow 
Miville-Deschenes et al. I d2008h model 4, using the same syn- 
chrotron spectral index map, and the synchrotron polarised tem- 
plate at 23 GHz measured by WMAP. This allows the definition 
of a pixel-dependent geometric depolarisation factor g(%), com- 
puted as the ratio between the polarisation expected theoretically 
from Eq.|2] and the polarisation actually observed. This depolar- 
isation, assumed to be due to varying orientations of the galactic 
magnetic field along the line of sight, is used also for modelling 
polarised dust emission (see below). 

As an additional refinement, we also investigate the impact 
of a slightly modified frequency dependence with a running 
spectral index in Sect. For this purpose, the synchrotron emis- 
sion Stokes parameters (Sf(g) for X e {Q, U}), at frequency v 
and in direction £ on the sky, will be modelled instead as: 



(3 



&(f)+C(f)log(v/v,) 



(3) 



where is the WMAP mea surement at v = 23GHz j3 s th e 

synchrotron spectral index map dMiville-Deschenes et al.l2 008). 
and C(£) a synthetic template of the curvature of the synchrotron 
spectral index. 

The reconstructed B-modes map of the synchrotron- 
dominated sky emission at 30 GHz is shown in Fig. [2] 

2.2. Dust 

The thermal emission from heated dust grains is the dominant 
galactic signal at frequencies higher than 100 GHz (Fig. [TJ. 
Polarisation of starlight by dust grains indicates partial align- 
ment of e longat ed grains with the galactic magnetic field (see 
iLazarianl d2007b for a review of possible alignment mecha- 
nisms). Partial alignment of grains should also result in polar- 
isation of the far infrared dust emission. 

Contributions from a wide range of grain sizes and composi- 
tions are required to expl ain the infrared spectrum of dust emis - 
sion from 3 to 1000 yum dDesert et alJll990tlLi & Drainell200lb . 
At long wavelengths of interest for CMB observations (above 
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100 fim), the emission from big grains, at equilibrium with the 
interstellar radiation field, should dominate. 

2.2.1. Dust thermal emission law 

There is no single theoretical emission law for dust, which is 
composed of many different populations of particles of mat- 
ter. On average, an emission law can be fit to observational 
data. In the frequency range of interest for CMB observations, 
Fink beiner et al.l (fl999) have shown that the dust emission in in- 
tensity is well modelled by emission from a two components 
mixture of silicate and carbon grains. For both components, the 
thermal emission spectrum is modelled as a modified grey-body 
emission, D v ~ B v (T)v a , with different emissivity spectral index 
a and different equilibrium temperature T. 

2.2.2. Dust polarisation 

So far, dust polarisation measurements have been mostly con- 
centrated on specific regions of emission, with the exceptio n 
of the Archeops balloon-borne experiment (Ben oit et al.l l2004). 
which has mapped the emission at 353 GHz on a significant part 
of the sky, showing a polarisation fraction around 4-5% and up 
to 10% in some clouds. This is in rough agreem ent with what 
could be expected from pola r isation of starlight (iFosalba et al.l 
l2002t iDraine & Fralssell2008l) . iMacellari et aTTd2008l) show that 
dust fractional polarisation in WMAP5 data depends on both fre- 
quency and latitude, but is typically about 3% and anyway below 

7% 

IDraine & Fraissd (|2008) have shown that for particular mix- 
tures of dust grains, the intrinsic polarisation of the dust emis- 
sion could vary significantly with frequency in the 100-800 GHz 
range. Geometrical depolarisation caused by integration along 
the line of sight also lowers the observed polarisation fraction. 

2.2.3. Our model of dust 

To summarise, dust produces polarised light depending on grains 
shape, size, composition, temperature and environment. The po- 
larised light is then observed after integration along a line of 
sight. Hence, the observed polarisation fraction of dust depends 
on its three-dimensional distribution, and of the geometry of the 
galactic magnetic field. This produces a complex pattern which 
is likely to be only partially coherent from one channel to an- 
other. 

Making use of the available data, the PSM models polarised 
thermal dust emission by extrapolating dust intensity to polar- 
isation intensity assuming an intrinsic polarisation fraction fa 
constant across frequencies. This value is set to fa = 0.12 
to be consistent w ith maximum values observed by Archeops 
(Benoit et al. 2004) and is in good agreement with the WMAP 
94 GHz measurement. The d ust intensity (pi), tr aced by the 
template m ap at 100 ^m fromlSchlegel et al.l dl998l) . is extrapo- 
lated using iFinkbeiner et al.l (1 19991 model #7) to frequencies of 
interest. The stokes Q and U parameters (respectively D Q and 
D u ) are then obtained as: 

= fd gig) D T v (g) cos(2y(£)) (4) 
D V v ig) =f d g(g)D T v ^sm(2y{g)) (5) 

The geometric 'depolarisation' factor g is a modified version of 
the synchrotron depolarisation factor (computed from WMAP 
measurements). Modifications account for differences of spa- 
tial distribution between dust grains and energetic electrons, 




-100.0 . . 100,0 uK 

Fig. 2. B-modes of the galactic foreground maps (synchrotron + 
dust) as simulated using vl.6.4 of the PSM. Top: synchrotron- 
dominated emission at 30 GHz, Bottom: dust-dominated emis- 
sion at 340 GHz. In spite of the fact that the direction of po- 
larisation of both processes is determined by the same galactic 
magnetic field, differences in the 3-D distributions and in the de- 
polarisation factors result in quite different B-mode polarisation 
patterns. 



and are computed using the m agnetic field model presented in 
Miville-Deschenes et al. (2008). The polarisation angle y is ob- 
tained from the magnetic field model on large scales and from 
synchrotron measurements in WMAP on scales smaller than 5 
degrees. Figure [2] shows the B-modes of dust at 340 GHz using 
this model. 

2.2.4. Anomalous dust 

If the anomalous dust emission, which may account for 
a significant part of the int e nsity emission in the range 
10-30 GHz dFinkbeineri l2004t Ide Oliveira-Costa et all 12004 
Miville-Deschenes et al.l 120081) . can be interpreted as s pin- 
ning dust grains emission (IDraine & Lazarianl 1 1998b, it 
shoul d be slightly polarised under 35 GHz (Battistell i et al.l 
2006 ), and only marginally polarised at higher frequencies 
(lLazarian & Finkbeiner! 2003). For this reason, it is neglected 
(and not modelled) here. However, we should keep in mind 
that there exist other possible emission processes for dust, like 
the magneto-dipole mechanism, which can produce highly po- 
larised radiation, and could thus contribute significantly to dust 
polarisation at low frequencie s, even if sub-dominant in intensity 
(Lazarian & Fink beiner 2003|). 

2.3. Other processes 

The left panel in Fig.Q]presents the respective contribution from 
the various foregrounds as predicted by the PSM at 100 GHz. 
Synchrotron and dust polarised emission, being by far the 
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strongest contaminants on large scales, are expected to be the 
main foregrounds for the measurement of primordial B-modes. 
In this work, we thus mainly focus on the separation from these 
two diffuse contaminants. However, other processes yielding po- 
larised signals at levels comparable with either the signal of in- 
terest, or with the sensitivity of the instrument used for B-mode 
observation, have to be taken into account. 

2.3.1. Free-free 

Free-free emission is assumed unpolarised to first order (the 
emission process is not intrinsically linearly polarised), even if, 
in principle, low level polarisation by Compton scattering could 
exist at the edge of dense io nised regions. In WMAP data anal- 
ysis, [Macellari et al] d2008l) find an upper limit of 1 % for free- 
free polarisation. At this level, free-free would have to be taken 
into account for measuring CMB B-modes for low values of r. 
As this is just an upper limit however, no polarised free-free is 
considered for the present work. 

2.3.2. Extra-galactic sources 

Polarised emission from extra-g alactic sources is e xpected to 
be faint below the degree scale. Tuccieta D (120051) . however, 
estimate that radio sources become the major contaminant af- 
ter subtraction of the galactic foregrounds. It is, hence, an im- 
portant foreground at high galactic latitudes. In addition, the 
point source contribution involves a wide range of emission 
processes and superposition of emissions from several sources, 
which makes this foreground poorly coherent across frequen- 
cies, and hence difficult to subtract using methods relying on the 
extrapolation of template emission maps. 

The Planck Sky Model provides estimates of the point 
source polarised e mission. Source cou nts are in agreement with 
the prediction of Ide Zotti et al.1 d2005l) . and with WMAP data. 
For radio-sources, the degree of polarisation for each source 
is randomly draw n from the observed distribution at 20 GHz 
dRicci et al.l2004l) . For infrared sources, a distribution with mean 
polarisation degree of 0.01 is assumed. For both populations, po- 
larisation angles are uniformly drawn in [0 - 2n[. The emission 
of a number of known galactic point sources is also included in 
PSM simulations. 

2.3.3. Lensing 

The last main contaminant to the primordial B-mode signal is 
lensing-induced B-type polarisation, the level of which should 
be of the same order as that of point sources (left panel of 
Fig. [TJ. For the present work, no sophisticated lensing clean- 
ing method is used. Lensing effects are modelled and taken into 
account only at the power spectrum level and computed using 
the CAM B software packaged based itself on the CMBFA ST 
software dZaldarriaga et al.lll998tlZaldarriaga & Seliakll2000l) . 

2.3.4. Polarised Sunyaev-Zel'dovich effect 

The polarised Sunyaev ZeFdovich effect jSazonov & Sunvaevl 
1999; I Audit & Simmons! 1 1 9991 ISeto & Pierpaolil 12005b . is ex- 
pected to be very sub-dominant and is neglected here. 



1 http://camb.info 



2.4. Uncertainties on the foreground model 

Due to the relative lack of experimental constraints from ob- 
servation at millimetre wavelengths, uncertainties on the fore- 
ground model are large. The situation will not drastically im- 
prove before the Planck mission provides new observations of 
polarised foregrounds. It is thus very important to evaluate, at 
least qualitatively, the impact of such uncertainties on compo- 
nent separation errors for B-mode measurements. 

We may distinguish two types of uncertainties, which impact 
differently the separation of CMB from foregrounds. One con- 
cerns the level of foreground emission, the other its complexity. 

Quite reliable constraints on the emission level of polarised 
synchrotron at 23 GHz are available with the WMAP measure- 
ment, up to the few degrees scale. Extrapolation to other fre- 
quencies and smaller angular scales may be somewhat insecure, 
but uncertainties take place where this emission becomes weak 
and sub-dominant. The situation is worse for the polarised dust 
emission, which is only weakly constrained from WMAP and 
Archeops at 94 and 353 GHz. The overall level of polarisation is 
constrained only in the galactic plane, and its angular spectrum 
is only roughl y estimated. In addition , variations of the polarisa- 
tion fraction dDraine & Fraissdl2008l) may introduce significant 
deviations to the frequency scaling of dust B-modes. 

Several processes make the spectral indexes of dust and syn- 
chrotron vary both in space and frequency. Some of this com- 
plexity is included in our baseline model, but some aspects, like 
the dependence of the dust polarisation fraction with frequency 
and the steepening of the synchrotron spectral index, remain 
poorly known and are not modelled in our main set of simula- 
tions. In addition, uncharacterised emission processes have been 
neglected. This is the case for anomalous dust, or polarisation 
of the free-free emission through Compton scattering. If such 
additional processes for polarised emission exist, even at a low 
level, they would decrease the coherence of galactic foreground 
emission between frequency channels, and hence our ability to 
predict the emission in one channel knowing it in the others - a 
point of much importance for any component separation method 
based on the combination of multi-frequency observations. 

The component separation as performed in this paper, hence, 
is obviously sensitive to these hypotheses. We will dedicate a 
part of the discussion to assess the impact of such modelling 
errors on our conclusions. 

3. Estimating / with contaminants 

Let us now turn to a presentation of the component separation 
(and parameter estimation) method used to derive forecasts on 
the tensor to scalar ratio measurements. 

Note that in principle, the best analysis of CMB observations 
should simultaneously exploit measu rements of all fields (T, 
E, and B), as investigated already by Au mont & Marias-P erez 
d2007l) . Their work, however, addresses an idealised problem. 
For component separation of temperature and polarisation to- 
gether, the best approach is likely to depend on the detailed prop- 
erties of the foregrounds (in particular on any differences, even 
small, between foreground emissions laws in temperature and 
in polarisation) and of the instrument (in particular noise cor- 
relations, and instrumental systematics). None of this is avail- 
able for the present study. For this reason, we perform compo- 
nent separation in B-mode maps only. Additional issues such 
as disentangling E from B in cases of partial sky coverage for 
instance, or in presence of instrumental systematic effects, are 
not investigated here either. Relevant work can be found in 



6 



Betoule et al.: T/S measurements in presence of foregrounds 



iKapkn &Delabrouilld j2002h : IChaUinor et all d2003b : iHu et al.l 
d2003l) : lRosset et alj (120071) . 

For low values of tensor fluctuations, the constraint on r 
is expected to come primarily from the B-mode polarisation. 
B-modes indeed are not affected by the cosmic variance of 
the scalar perturbations, contrarily to E-modes and temperature 
anisotropies. In return, B-mode signal would be low and should 
bring little constraint on cosmological parameters other than r 
(and, possibly, the tensor spectral index «,, although this addi- 
tional parameter is not considered here). Decoupling the esti- 
mation of r (from B-modes only) from the estimation of other 
cosmological parameters (from temperature anisotropies, from 
E-modes, and from additional cosmological probes) thus be- 
comes a reasonable hypothesis for small values of r. As we 
are primarily interested in accurate handling of the foreground 
emission, we will make the assumption that all cosmological 
parameters but r are perfectly known. Further investigation of 
the coupling between c osmological param eters can be found in 
IColombo et all d2008l) : IVerde et al.l (120061) . and this question is 
discussed a bit further in Sect. 



3.1.2. Multi-map estimation 

Alternatively, we may consider observations in F frequency 
bands, and form the Fx 1 vector of data x(£), assuming that 
each frequency is contaminated by jc cont . This term includes all 
contaminations (foregrounds, noise, etc.). In the harmonic do- 
main, denoting A cm b the emission law of the CMB (the unit vec- 
tor when working in thermodynamic units): 



„ —A ^, cm b i „cm« 



(9) 



We then consider the F x F spectral covariance matrix R e con- 
taining auto and cross-spectra. The CMB signal being uncorre- 
cted with the contaminants, one has: 



with the CMB contribution modelled as 



R C ™V) = rS c A cmh Al 



(10) 



(11) 



3.1. Simplified approaches 
3.1.1. Single noisy map 

The first obstacle driving the performance of an experiment be- 
ing the instrumental noise, it is interesting to recall the limit on 
r achievable in absence of foreground contamination in the ob- 
servations. 

We thus consider first a single frequency observation of the 
CMB, contaminated by a noise term n: 

x(f) = x™ b ({) + „(0 (6) 

where £ denotes the direction in the sky. Assuming that n is un- 
corrected with the CMB, the power spectra of the map reads: 

C( = rS ( + N t 

where St is the shape of the CMB power-spectrum (as set by 
other cosmological parameters), and Nt the power of the noise 
contamination. Neglecting mode to mode mixing effects from 
a mask (if any), or in general from incomplete sky coverage, 
and assuming that n can be modelled as a Gaussian process, the 
log-likelihood function for the measured angular power spec- 
trum reads: 



~21n£= Y(2e + l)/ 8ky ln(^) + ^ 
i Y \Cel w 



+ const. 



(7) 



The smallest achievable variance cr 2 r in estimating r is the 
inverse of the Fisher information I = -Ef which takes 

the form: 



cr. 



-2 



2t+ 1 



/sky 



rSc + Nt 



(8) 



For a detector (or a set of detectors at the same frequency) 
of noise equivalent temperature s (in fiK -\fs), and a mission 
duration of t s seconds, the detector noise power spectrum is 
Nt = ^If- /xK 2 , with Be denoting the beam transfer function of 
the detector. 

A similar approach to estimating o> is used in IVerde et ail 
(2006) where a single 'cleaned' map is considered. This map is 
obtained by optimal combination of the detectors with respect to 
the noise and cleaned from foregrounds up to a certain level of 
residuals, which are accounted for as an extra Gaussian noise. 



and all contaminations contributing a term Nf to be discussed 
later. The dagger (f) denotes the conjugate transpose for com- 
plex vectors and matrices, and the transpose for real matrices (as 

-^cmb)- 

In the approximation that contaminants are Gaussian (and, 
here, stationary) but correlated, all the relevant information 
about the CMB is preserved by combining all the channels into 
a single filtered map. In the harmonic domain, the filtering oper- 
ation reads: 

W ( a c „ 



aim 



a% b +W t a 



with 



^cmb ^ 



cont 



A Lb N ^' A cmb 



(12) 



We are back to the case of a single map contaminated by a 
characterised noise of spectrum: 



Nt = E\Wta?«f = (Al mb Nt- 1 A c J ) ~ 



(13) 



If the residual WVa?L nt is modelled as Gaussian, the single- 
map likelihood © can be used 

The same filter is used bv lAmblard et al.l (120071) . Assuming 
that the foreground contribution is perfectly known, the contam- 
inant terms Nf can be modelled as Nf = R" olse + R^ E . This ap- 
proach thus permits to derive the actual level of contamination 
of the map in presence of known foregrounds, i.e. assuming that 
the covariance matrix of the foregrounds is known. 



3.2. Estimating r in presence of unknown foregrounds with 
SMICA 



The two simplified approaches of sections 13.1.11 and 13.1.21 offer 
a way to estimate the impact of foregrounds in a given mission, 
by comparing the sensitivity on r obtained in absence of fore- 
grounds (from Eq.[8]when Nt contains instrumental noise only), 
and the sensitivity achievable with known foregrounds (when Nt 
contains the contribution of residual contaminants as well, as 
obtained from Eq. [13] assuming that the foreground correlation 
matrix is known). 

A key issue, however, is that the solution and the error bar 
require the covariance matrix of foregrounds and noise to be 
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knownQ Whereas the instrumental noise can be estimated accu- 
rately, assuming prior knowledge of the covariance of the fore- 
grounds to the required precision is optimistic. 

To deal with unknown foregrounds, we thus follow 
a different route which considers a multi-map likelihood 
dDelabrouille et all 120031) . If all processes are modelled as 
Gaussian isotropic, then standard computations yield: 

-21n£ = + l)/sky^(R f , R f ) + est (14) 



where FL is the sample estimate of FL : 



It 



--E 

+ 1 f s kv A 



a Um a lm 



(15) 



m=-{ 



and where K (•, •) is a measure of mismatch between two positive 
matrices given by: 

K(R, r) = - [trace(Fr'R) - log det(R _1 R) - f] (16) 

Expression dT4l > is nothing but the multi-map extension of d7j. 

If Nf is known and fixed, then the likelihood (Eq. IT4T > de- 
pends only on the CMB angular spectrum and can be shown to 
be equal (up to a constant) to expression [7] with Q = rS t and 
N{ given by Eq. Qj] Thus this approach encompasses both the 
single map and filtered map approaches. 

Unknown foreground contribution can be modelled as the 
mixed contribution of D correlated sources: 

R[ g = AZ f A + (17) 

where A is a F x D mixing matrix and If is the D x D spec- 
tral covariance matrix of the sources. The model of the spectral 
covariance matrix of the observations is then: 

R, = rS ( A cmb Al mb + AI,A + + R™ ise 

We then maximise the likelihood ( fT4b of the model with respect 
to r, A and If. 

We note that the foreground parameterisation in Eq[l7]is re- 
dundant, as an invertible matrix can be exchanged between A 
and I, without modifying the actual value of R fg . The physical 
meaning of this is that the various foregrounds are not identi- 
fied and extracted individually, only their mixed contribution is 
characterised. 

If we are interested in disentangling the foregrounds as 
well, e.g. to separate synchrotron emission from dust emission, 
this degeneracy can be lifted by making use of prior informa- 
tion to constrain, for example, the mixing matrix. Our multi- 
dimensional model offers, however, greater flexibility. Its main 
advantage is that no assumption is made about the foreground 
physics. It is not specifically tailored to perfectly match the 
model used in the simulation. Because of this, it is generic 
enough to absorb variations in the properties of the foregrounds, 
as will be seen later-on, but specific enough to preserve iden- 
tifiability in the separation of CMB from foreground emission. 
A more complete discussion of the Smica metho d with flexible 
components can be found in lCardoso et alj J20Q8). 

A couple last details on Smica and its practical implemen- 
tation are of interest here. For numerical purposes, we actu- 
ally divide the whole I range into Q frequency bins T> q = 



2 The actual knowledge of the contaminant term is not strictly re- 
quired to build the filter. It is required, however, to derive the contami- 
nation level of the filtered map. 



K™ n , ■ ■ ■ ,^™ ax ), and form the binned versions of the empirical 
and true cross power-spectra: 



_ 1 e 

H « = - Z Z ai ' ma l 

i eeD„m=-t 

R ? =^Z^ +1 ^ 



(18) 



CeD a 



where w q is the number of modes in D q . It is appropriate to 
select the domains so that we can reasonably assume for each 
£ € D q , R f ~ R„. This means that spectral bins should be small 
enough to capture the variations of the power spectra. In prac- 
tice results are not too sensitive to the choice of the spectral bin 
widths. Widths between 5 and 10 multipoles constitute a good 
tradeoff. 

Finally, we compute the Fisher information matrix \jj(0) de- 
riving from the maximised likelihood ( TBI for the parameter set 
= (r,A,Zi,---,ZQ): 



km 



w„ trace 



(dR q (ff),dR q (ff), 



\ d6i 



(19) 



The lowest achievable variance of the r estimate is obtained as 
the entry of the inverse of the FIM corresponding to the param- 
eter r. 



-l 



(20) 



4. Predicted results for various experimental 
designs 

We now turn to the numerical investigation of the impact of 
galactic foregrounds on the measurements of r with the follow- 
ing experimental designs: 

- The Planck space mission, due for launch early 2009, which, 
although not originally planned for B-mode physics, could 
provide a first detection if the tensor to scalar ratio r is around 
0.1. 

- Various versions of the EPIC space mission, either low cost 
and low resolution (EPIC-LC), or more ambitious versions 
(EPIC-CS and EPIC-2m). 

- An ambitious (fictitious) ground-based experiment, based on 
the extrapolation of an existing design (the Clover experi- 
ment). 

- An alternative space mission, with sensitivity performances 
similar to the EPIC-CS space mission, but mapping only a 
small (and clean) patch of the sky, and referred as the 'deep 
field mission'. 

The characteristics of these instruments are summed-up in ta- 
ble [T] and Fig. [3] illustrates their noise angular power spectra in 
polarisation. 

4.1. Pipeline 

For each of these experiments, we set up one or more simula- 
tion and analysis pipelines, which include, for each of them, the 
following main steps: 

- Simulation of the sky emission for a given value of r and a 
given foreground model, at the central frequencies and the 
resolution of the experiment. 
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Planck 

EPIC-LC 

EPIC-CS 

EPIC-2m 

Ground-Based 



Planck> 




Fig. 3. Noise spectra of various experimental designs compared 
to B-modes levels for r = 0.1, 0.01 and 0.001. When computing 
the equivalent multipole noise level for an experiment, we as- 
sume that only the central frequency channels contribute to the 
CMB measurement and that external channels are dedicated to 
foreground characterisation. 




Fig. 4. Analysis mask for EPIC B maps, smoothed with a l c 
apodisation window. 



Simulation of the experimental noise, assumed to be white, 
Gaussian and stationary. 

Computation, for each of the resulting maps, of the coeffi- 
cients of the spherical harmonic expansion of the B-modes 

a B r 

tm 

Synthesis from those coefficients of maps of B-type signal 
only. 

For each experiment, a mask based on the B-modes level 
of the foregrounds is built to blank out the brightest fea- 
tures of the galactic emission (see Fig.|4]i. This mask is built 
with smooth edges to reduce mode-mixing in the pseudo- 
spectrum. 

Statistics described in Equation[18]are built from the masked 
B maps. 

The free parameters of the model described in Sect. 13.21 are 
adjusted to fit these statistics. The shape of the CMB pseudo- 
spectrum that enters in the mod el, is computed u sing the 
mode-mixing matrix of the mask dHivon et al.ll2002l) . 
Error bars are derived from the Fisher information matrix of 
the model. 



Some tuning of the pipeline is necessary for satisfactory 
foreground separation. The three main free parameters are the 
multipole range [Am n ,AnaxL the dimension D of the foreground 
component, and (for all-sky experiments) the size of the mask. 

In practice we choose i m [ n according to the sky coverage and 
Anax according to the beam and the sensitivity. The value of D 
is selected by iterative increments until the goodness of fit (as 
measured from the Smica criterion on the data themselves, with- 
out knowledge of the input CMB and foregrounds) reaches its 
expectation. The mask is chosen in accordance to maximise the 
sky coverage for the picked value of D (see appendix lAl for fur- 
ther discussion of the procedure). 

For each experimental design and fiducial value of r we com- 
pute three kinds of error estimates which are recalled in Table [2] 

Knowing the noise level and resolution of the instrument, we 



first derive from Eq. [8] the error cr" 



noise-only 



set by the instrument 



sensitivity assuming no foreground contamination in the covered 
part of the sky. The global noise level of the instrument is given 

by N[ = (A^N^'Acmb) , where the only contribution to Nf 

/ 4;r 4 \ 

comes from the instrumental noise: H( = R" olse = diag I gf p t I- 

T ,« i ,i known-foreground 

In the same way, we also compute the error o> 
that would be obtained if foreground contribution R fg to the co- 
variance of the observations was perfectly known, using Nf = 
R° oise + R^ s . Here we assume that R fg = R fg where R fg is the 
sample estimate of R fg computed from the simulated foreground 
maps. 

Finally, we compute the error cr^ MICA given by the Fisher 
information matrix of the model (Eq.l20l. 

In each case, we also decompose the FIM in the contribution 
from large scale modes (i < 20) and the contribution from small 
scales (£ > 20) to give indications of the relative importance of 
the bump (due to reionisation) and the peak (at higher €) in the 
constraint of r. 

We may notice that in some favourable cases (at low I, where 
the foregrounds dominate), the error estimate given by Smica can 
be slightly more optimistic than the estimate obtained using the 
actual empirical value of the correlation matrix R fg . This reflects 
the fact that our modelling hypothesis, which imposes to R tg to 
be of rank smaller than D, is not perfectly verified in practice 
(see Appendix lAl for further discussion of this hypothesis). The 
(small) difference (an error on the estimation of cr r when fore- 
grounds are approximated by our model) has negligible impact 
on the conclusions of this work. 

4.2. Planck 

The Planck space mission will be the first all-sky experiment to 
give sensitive measurements of the polarised sky in seven bands 
between 30 and 353 GHz. The noise level of this experiment be- 
ing somewhat too high for precise measurement of low values of 
r, we run our pipeline for r — 0.1 and 0.3. We predict a possible 
3-sigma measurement for r — 0.1 using Smica (first lines in ta- 
ble|2]). A comparison of the errors obtained from Smica, with the 
prediction in absence of foreground contamination, and with per- 
fectly known foreground contribution, indicates that the error is 
dominated by cosmic variance and noise, foregrounds contribut- 
ing to a degradation of the error of ~ 30% and uncertainties on 
foregrounds for another increase around 30% (for r = 0.1). 

Fig.|3]hints that a good strategy to detect primordial B-modes 
with Planck consists in detecting the reionisation bump below 
I = 10, which requires the largest possible sky coverage. Even 
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Experiment 


frequency 


beam FWHM 


NET 


Tobs 


sky coverage 




(GHz) 


(') 




(yr) 


(fsky) 


PLANCK 


30, 44, 70 


33, 24, 14 


96, 97, 97 


1.2 


1 


100, 143, 217, 353 


10, 7.1, 5, 5 


At O 1 C 1 1 C A 

41, 31, 51, 154 






EPIC-LC 


30, 40, 60 


155, 116, 77 


28, 9.6, 5.3 


2 


1 


dpi lie inn onn 

90, 135, 200, 300 


52, 34, 23, 16 


O T 'IT TO 

2.3, 2.2, 2.3, 3.8 






EPIC-CS 


30, 45, 70, 100 
150, 220, 340, 500 


15.5, 10.3,6.6, 4.6 
3.1,2.1, 1.4, 0.9 


19, 8, 4.2,3.2 
3.1, 5.2, 25, 210 


4 


1 


EPIC-2m 


30, 45, 70, 100 
150, 220, 340, 500(,800) 


26, 17, 11, 8 
5,3.5, 2.3, 1.5(, 0.9) 


18, 7.6,3.9,3.0 
2.8, 4.4, 20, 180(, 28k) 


4 


1 


Ground-Based 


97, 150, 225 


7.5,5.5,5.5 


12, 18, 48 


0.8 


0.01 


Deep field 


30, 45, 70, 100 


15.5, 10.3,6.6, 4.6 


19, 8, 4.2,3.2 


4 


0.01 


150, 220, 340, 500 


3.1,2.1, 1.4, 0.9 


3.1, 5.2, 25, 210 







Table 2. Error prediction for various experimental designs and fiducial r values. Error bars from the columns noise-only and known 
foregrounds are derived from Eq.© assuming Nf = R nolse and H( = R nolse + R f s respectively. Error bars from the Smica column are 
obtained by the inversion of the FIM computed from the Smica model at the point of convergence of the algorithm as in Eq.(l20ll. 
In all cases, large scale (I < 20) and small scale {£ > 20) error bars are computed by decomposing the Fisher information between 
contribution from low and high multipoles. This allows for an estimation of respective contribution from the bump and the peak 
to the measurement. The r est column gives the estimated value at the convergence point in Smica. Detections at more than Act are 
bold-faced. 



case 


r 


(T T /r 


noise-only 

of 20 /r of 20 /r 


known foregrounds 
a-r/r of 20 /r a- e r >20 /r 


<r r /r 


Smica 

of 20 /r 


of 20 /r 


r est 


Anin 


^max 


/sky 


D 3 


PLANCK 


0.3 


0.075 


0.17 


0.084 


0.1 


0.2 


0.12 


0.15 


0.22 


0.2 


0.26 




130 


0.95 




0.1 


0.17 


0.25 


0.22 


0.23 


0.34 


0.32 


0.29 


0.34 


0.55 


0.086 


2- 


3 


EPIC-LC 


0.01 


0.019 


0.084 


0.019 


0.05 


0.18 


0.053 


0.079 


0.18 


0.1 


0.0098 


2- 


130 


0.86 


4 


0.001 


0.059 


0.15 


0.064 


0.27 


0.4 


0.38 


0.37 


0.43 


0.82 


0.00088 






EPIC-2m 


0.01 
0.001 


0.016 
0.051 


0.083 
0.14 


0.016 
0.055 


0.027 
0.14 


0.12 
0.25 


0.027 
0.16 


0.032 
0.16 


0.11 
0.24 


0.036 
0.24 


0.0096 
0.001 


2- 


300 


0.87 


4 


EPIC-CS 


0.01 


0.017 


0.084 


0.017 


0.029 


0.12 


0.03 


0.036 


0.11 


0.041 


0.0096 


2- 


300 


0.87 


4 


0.001 


0.058 


0.15 


0.063 


0.15 


0.27 


0.19 


0.18 


0.26 


0.29 


0.00098 






Ground-based 


0.1 
0.01 


0.083 
0.18 






0.15 
0.8 






0.24 
1.6 






0.11 

0.018 


50 


-300 


0.01 


2 


Grnd-based+Planck 


0.01 


0.18 


0.51 






0.69 






0.0065 


50 


-300 


0.01 


2 


Deep field mission 


0.001 


0.082 






0.1 






0.13 






0.00092 


50 


-300 


0.01 
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at high latitude, a model using D — 2 fails to fit the galactic 
emission, especially on large scales where the galactic signal is 
above the noise. Setting D — 3, however, gives a satisfactory 
fit (as measured by the mismatch criterion) on 95 percent of the 
sky. It is therefore our choice for Planck. 

We also note that a significant part of the information is com- 
ing from the reionisation bump (£ < 20). The relative importance 
of the bump increases for decreasing value of r, as a consequence 
of the cosmic variance reduction. For a signal-to-noise ratio cor- 
responding roughly to the detection limit (r = 0.1), the stronger 
constraint is given by the bump (Appendix [B] gives further il- 
lustration of the relative contribution of each multipole). This 
has two direct consequences: the result is sensitive to the actual 
value of the reionis ation optical depth and to reio nisation history 
(as investigated bv lColom bo & Pierpaoli (2008)), and the actual 
capability of Planck to measure r will depend on the level (and 
the knowledge of) instrumental systematics on large scales. 

Note that this numerical experiment estimates how well 
Planck can measure r in presence of foregrounds from B-modes 
only. 

4.3. EPIC 

We perform a similar analysi s for three possible designs of the 
EPIC probe dBock et al.ll2008h . EPIC-LC and EPIC-CS corre- 



spond respectively to the low cost and comprehensive solutions. 
EPIC-2m is an alternate design which contains one extra high- 
frequency channel (not considered in this study) dedicated to ad- 
ditional scientific purposes besides CMB polarisation. We con- 
sider two values of r, 0.01 and 0.001. For all these three experi- 
ments, the analysis requires D = 4 for a reasonable fit, which is 
obtained using about 87% of the sky. 

The two high resolution experiments provide measurements 
of r = 10~ 3 with a precision better than five sigma. For the lower 
values of r, the error is dominated by foregrounds and their pres- 
ence degrades the sensitivity by a factor of 3, as witnessed by 
the difference between cr ™ lse onl >' an( j cr™ lca . However, while the 
difference between the noise-only and the Smica result is a fac- 
tor 4-6 for EPIC-LC, it is only a factor about 2-3 for EPIC-CS 
and EPIC-2m. Increased instrumental performance (in terms of 
frequency channels and resolution) thus also allows for better 
subtraction of foreground contamination. 

For all experiments considered, the constraining power 
moves from small scales to larger scale when r decreases down 
to the detection limit of the instrument. In all cases, no informa- 
tion for the CMB is coming from t > 150. Higher multipoles, 
however, are still giving constraints on the foreground param- 
eters, effectively improving the component separation also on 
large scales. 
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4.4. Small area experiments 
AAA. Ground-based 

A different observation strategy for the measurement of B-modes 
is adopted for ground-based experiments that cannot benefit 
from the frequency and sky coverage of a space mission. Such 
experiments target the detection of the first peak around I = 100, 
by observing a small but clean area (typically 1000 square- 
degrees) in few frequency bands (2 or 3). 

The test case we propo se here is inspired from the announced 
performances of Clover dNorth et al.ll2008h . The selected sky 
coverage is a 10 degree radius area centred on Ion = 351°, 
lat = -56° in galactic coordinates. The region has been retained 
by the Clover team as a tradeoff between several issues includ- 
ing, in particular, foreground and atmospheric contamination. 
According to our polarised galactic foreground model, this also 
correspond to a reasonably clean part of the sky (within 30% of 
the cleanest). 

The most interesting conclusion is that for r = 0.01, although 
the raw instrumental sensitivity (neglecting issues like E-B mix- 
ing due to partial sky coverage) would allow a more than five 
sigma detection, galactic foregrounds cannot be satisfactorily re- 
moved with the scheme adopted here. 

An interesting option would be to complement the measure- 
ment obtained from the ground, with additional data as that of 
Planck, and extract r in a joint analysis of the two data sets. 
To simply test this possibility here, we complement the ground 
data set with a simulation of the Planck measurements on the 
same area. This is equivalent to extend the frequency range of the 
ground experiment with less sensitive channels. We find a signif- 
icant improvement of the error-bar from 1 .6 ■ 10~ 2 to 0.69 ■ 10 -2 , 
showing that a joint analysis can lead to improved component 
separation. The degradation of sensitivity due to foreground re- 
mains however higher than for a fully sensitive space mission (as 
witnessed by the following section). This last result is slightly 
pessimistic as we do not make use of the full Planck data set but 
use it only to constrain foregrounds in the small patch. However 
considering the ratio of sensitivity between the two experiments, 
it is likely that there is little to gain by pushing the joint analysis 
further. 

4.4.2. Deep field space mission 

We may also question the usefulness of a full-sky observation 
strategy for space-missions, and consider the possibility to spend 
the whole observation time mapping deeper a small but clean 
region. 

We investigate this alternative using an hypothetical experi- 
ment sharing the sensitivity and frequency coverage of the EPIC- 
CS design, and the sky coverage of the ground-based experi- 
ment. Although the absence of strong foreground emission may 
permit a design with a reduced frequency coverage, we keep a 
design similar to EPIC-CS to allow comparisons. In addition, the 
relative failure of the ground-based design to disentangle fore- 
grounds indicates that the frequency coverage cannot be freely 
cut even when looking in the cleanest part of the sky. In the same 
way, to allow straightforward comparison with the ground-based 
case we stick to the same sky coverage, although in principle, 
without atmospheric constraints, slightly better sky areas could 
be selected. 

In spite of the increased cosmic variance due to the small 
sky coverage, the smaller foreground contribution allows our 
harmonic-based foreground separation with Smica to achieve 



better results with the 'deep field' mission than with the full sky 
experiment, when considering only diffuse galactic foreground. 
However, this conclusion doesn't hold if lensing is considered as 
will be seen in the following section. 

We may also notice that, despite the lower level of fore- 
grounds, the higher precision of the measurement requires the 
same model complexity (D = 4) as for the full sky experiment 
to obtain a good fit. 

We also recall that our processing pipeline does not ex- 
ploit the spatial variation of foreground intensity, and is, in this 
sense, suboptimal, in particular for all-sky experiments. Thus, 
the results presented for the full-sky experiment are bound to 
be slightly pessimistic which tempers further the results of this 
comparison between deep field and full sky mission. This is fur- 
ther discussed below. 

Finally, note that here we also neglect issues related to partial 
sky coverage that would be unavoidable in this scheme. 

4.5. Comparisons 

4.5.1 . Impact of foregrounds: the ideal case 

As a first step, the impact of foregrounds on the capability to 
measure r with a given experiment, if foreground covariances 
are known, is a measure of the adequacy of the experiment to 
deal with foreground contamination. Figures for this comparison 
are computed using equations [8] and [13] and are given in table [2] 
(first two sets of three columns). 

The comparison shows that for some experiments, cr r /r in 
the 'noise-only' and the 'known foregrounds' cases are very 
close. This is the case for Planck and for the deep field mis- 
sion. For these experiments, if the second order statistics of 
the foregrounds are known, galactic emission does not impact 
much the measurement. For other experiments, the 'known fore- 
grounds' case is considerably worse than the 'noise-only' case. 
This happens, in particular, for a ground based experiment when 
r = 0.01, and for EPIC-LC. 

If foreground contamination was Gaussian and stationary, 
and in absence of priors on the CMB power spectrum, the linear 
filter of equation 1121 would be the optimal filter for CMB re- 
construction. The difference between <j r in the 'noise-only' and 
the 'known foregrounds' cases would be a good measure of how 
much the foregrounds hinder the measurement of r with the ex- 
periment considered. A large difference would indicate that the 
experimental design (number of frequency channels and sensi- 
tivity in each of them) is inadequate for 'component separation' . 

However, since foregrounds are neither Gaussian nor sta- 
tionary, the linear filter of equation [12] is not optimal. Even if 
we restrict ourselves to linear solutions, the linear weights given 
to the various channels should obviously depend on the local 
properties of the foregrounds. Hence, nothing guarantees that we 
can not deal better with the foregrounds than using a linear fil- 
ter in harmonic space. Assuming that the covariance matrix of 
the foregrounds is known, the error in equation [8] with Nc from 
equation [T3l is a pessimistic bound on the error on r. The only 
conclusion that can be drawn is that the experiment does not al- 
low effective component separation with the implementation of 
a linear filter in harmonic space. There is, however, no guaran- 
tee either that an other approach to component separation would 
yield better results. 

Hence, the comparison of the noise-only and known fore- 
grounds cases shown here gives an upper limit of the impact of 
foregrounds, if they were known. 
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4.5.2. Effectiveness of the blind approach 

Even if in some cases the linear filter of equation Q~2] may not 
be fully optimal, it is for each mode £ the best linear com- 
bination of observations in a set of frequency channels, to re- 
ject maximally contamination from foregrounds and noise, and 
minimise the error on r. Other popular methods as decorrela- 
tion in direct space, as the so-called 'internal linear combina- 
tion' (ILC), and other linear combinations cannot do better, un- 
less they are implemented locally in both pixel and harmonic 
space simultaneously, us ing for instance spherical needlets as in 
Dela brouille et al.l |2008). Such localisation is not considered in 
the present work. 

Given this, the next question that arises is how well the spec- 
tral covariance of the foreground contamination can be actually 
constrained from the data, and how this uncertainty impact the 
measurement of r. The answer to this question is obtained by 
comparing the second and third sets of columns of table [2] 

In all cases, the difference between the results obtained as- 
suming perfect knowledge of the foreground residuals, and those 
obtained after the blind estimation of the foreground covariances 
with Smica, are within a factor of 2. For EPIC-2m and the deep 
field mission, the difference between the two is small, which 
means that Smica allows for component separation very effec- 
tively. For a ground based experiment with three frequency chan- 
nels, the difference is very significant, which means that the data 
does not allow a good blind component separation with Smica. 

Comparing column set 1 (noise-only) and 3 (blind approach 
with Smica) gives the overall impact of unknown galactic fore- 
grounds on the measurement of r from B-modes with the var- 
ious instruments considered. For Planck, EPIC-2m, or a deep 
field mission with 8 frequency channels, the final error bar on 
r is within a factor of 2 of what would be achievable without 
foregrounds. For EPIC-LC, or even worse for a ground-based 
experiment, foregrounds are likely to impact the outcome of the 
experiment quite significantly. For this reason, EPIC-2m and the 
deep field mission seem to offer better perspectives for measur- 
ing r in presence of foregrounds. 

4.5.3. Full sky or deep field 

The numerical investigations performed here allow -to some 
extent- to compare what can be achieved with our approach in 
two cases of sky observation strategies with the same instrument. 
For EPIC-CS, it has been assumed that the integration time is 
evenly spread on the entire sky, and that 87% of the sky is used 
to measure r. For the 'deep field' mission, 1 % of the sky only is 
observed with the same instrument, with much better sensitivity 
per pixel (by a factor of 10). 

Comparing cr r /r between the two in the noise-only case 
shows that the full sky mission should perform better (by a factor 
1 .4) if the impact of the foregrounds could be made to be negligi- 
ble. This is to be expected, as the cosmic or 'sample' variance of 
the measurement is smaller for larger sky coverage. After com- 
ponent separation however, the comparison is in favour of the 
deep field mission, which seems to perform better by a factor 
1.4 also. The present work, however, does not permit to con- 
clude on what is the best strategy for two reasons. First, this 
study concentrates on the impact of diffuse galactic foregrounds 
which are not expected to be the limiting issue of the deep field 
design. And secondly, in the case of a deep field, the properties 
of the (simulated) foreground emission are more homogeneous 
in the observed area, and thus the harmonic filter of equation[T2l 
is close to optimal everywhere. For the full sky mission, how- 



ever, the filter is obtained as a compromise minimising the over- 
all error £ by £, which is not likely to be the best everywhere on 
the sky. Further work on component separation, making use of a 
localised version of Smica, is needed to conclude on this issue. 
A preliminary version of Smica in wavelet space is described in 
Mou dden et alJd2004l) . but applications to CMB polarisation and 
full sky observations require specific developments. 

5. Discussion 

The results presented in the previous section have been obtained 
using a number of simplifying assumptions. 

First of all, only galactic foregrounds (synchrotron and dust) 
are considered. It has been assumed that other foregrounds (point 
sources, lensing) can be dealt with independently, and thus will 
not impact much the overall results. 

Second, it is quite clear that the results may depend on de- 
tails of the galactic emission, which might be more complex than 
what has been used in our simulations. 

Third, most of our conclusions depend on the accuracy of 
the determination of the error bars from the Fisher information 
matrix. This method, however, only provides an approximation, 
strictly valid only in the case of Gaussian processes and noise. 

Finally, the measurement of r as performed here assumes 
a perfect prediction (from other sources of information) of the 
shape of the BB spectrum. 

In this section, we discuss and quantify the impact of these 
assumptions, in order to assess the robustness of our conclu- 
sions. 

5.1. Small scale contamination 
5.1.1. Impact of lensing 

Limitations on tensor mode detection due to lensing have been 
widely investigated in the literature, and cleaning methods, 
based on the reconstruction of the lensed B-modes from es- 
timation of the len s potential and unlensed CMB E-modes , 
have been proposed (Knox & Song 2002; Hirata & Seliak 2003; 
iKesden et al] l2003b iLewis & Challinorll2006l) . However, limits 
on r achievable after such 'delensing' (if any) are typically sig- 
nificantly lower than limits derived in Sect. |H for which fore- 
grounds and noise dominate the error. 

In order to check whether the presence of lensing can signifi- 
cantly alter the detection limit, we proceed as follows: assuming 
no specific reconstruction of the lens potential, we include lens- 
ing effects in the simulation of the CMB (at the power spectrum 
level). The impact of this on the second order statistics of the 
CMB is an additional contribution to the CMB power spectrum. 
This extra term is taken into account on the CMB model used 
in Smica. For this, we de-bias the CMB Smica component from 
the (expectation value of) the lensing contribution to the power- 
spectrum. The cosmic variance of the lensed modes thus con- 
tributes as an extra 'noise' which lowers the sensitivity to the pri- 
mordial signal, and reduces the range of multipoles contributing 
significantly to the measurement. We run this lensing test case 
for the EPIC-CS and deep field mission. Table [3] shows a com- 
parison of the constraints obtained with and without lensing in 
the simulation for a fiducial value of r — 0.001. On large scales 
for EPIC-CS, lensing has negligible impact on the measurement 
of r (the difference between the two cases, actually in favour of 
the case with lensing, is not significant on one single run of the 
component separation). On small scales, the difference becomes 
significant. Overall, o~ r /r changes from 0.17 to 0.2, not a very 



12 



Betoule et al.: T/S measurements in presence of foregrounds 



significant degradation of the measurement: lensing produces a 
15% increase in the overall error estimate, the small scale error 
(for I > 20) being most impacted. For the small coverage mis- 
sion, however, the large cosmic variance of the lensing modes 
considerably hinder the detection. 







no lensing 




lensing 


Experiment 


<r r /r 


of^/r cr r r >20 /r 


o>/r 


<r^ 20 /r o- c > 20 /r 


EPIC-CS 


0.17 


0.25 0.28 


0.2 


0.24 0.36 


Deep field 


0.13 




1.1 





Table 3. Comparison of the constraints on r with and without 
lensing (here r = 0.001). 



Thus, at this level of r, if the reionisation bump is satisfac- 
torily measured, the difference is perceptible but not very sig- 
nificant. Hence, lensing is not the major source of error for a 
full-sky experiment measuring r. It becomes however a poten- 
tial problem for a small coverage experiment targeting the mea- 
surement of the recombination bump. Such a strategy would thus 
require efficient 'delensing'. Indications that 'delensing' can be 
performed even in presence of foregrounds in the c ase of a low 
noise a nd high resolution experiment can be found in Smit h et all 
(2008). However, a complete investigation of this case, account- 
ing for all the complexity (diffuse foregrounds, point sources, 
lensing, modes-mixing effects), would be needed to conclude on 
the validity of a deep-field strategy. 

5.1 .2. Impact of extra-galactic sources 

Although largely sub-dominant on scales larger than 1 degree, 
extra-galactic sources, in particular radio-sources, a re expected 
to be t h e worst contaminant on sm all scales (see e.g. lTucci et alj 
(2004); lPierpaoli & Pernal ( 120041) ). 

Obviously, the strongest point sources are known, or (for 
most of them) will be detected by Planck. Their polarisation 
can be measured either by the B-mode experiment itself, or by 
dedicated follow-up. We make the assumption that point sources 
brighter than 500 mJy in temperature (around 6000 sources) are 
detected, and that their polarised emission is subtracted from the 
polarisation observations. We stress that 500 mJy is a conserva- 
tive assumption as Planck is expected to have better detection 
thresholds. 

The present level of knowledge about point sources does not 
allow a very accurate modelling of the contribution to the power 
spectra of the remaining point sources (those not subtracted by 
the 500 mJy cut). For this reason we investigate their impact 
in two extreme cases: perfect modelling of their contribution to 
the power-spectra ('ideal' case), and no specific modelling at 
all ('no-model' case). Results of a Smica run for both assump- 
tions are compared to what is obtained in total absence of point 
sources ('no-ps' case), and are summarised in table|4] 

The bottom line of this investigation is that modelling prop- 
erly the point sources statistical contribution is necessary to mea- 
sure r = 0.001. An insufficient model results in a biased es- 
timator: for EPIC-CS the estimated r is two times larger than 
expected, with a difference incompatible with the error bar, in 
spite of an increased standard deviation (cry increased by +30% 
for r = 0.001). 

An ideal model restores the goodness of fit of the no-ps case 
and suppresses the bias of the estimator. Still, the presence of 




100 150 200 250 300 



Fig. 5. Goodness-of-fit for the three point sources cases. For the 
reference case 'no-ps', point sources have neither been including 
in the simulation, nor taken into account in the modelling. The 
mismatch criterion wander around its expectation value (hori- 
zontal dashed line). The 'no-model' case is a pessimistic situa- 
tion where no effort has been made to model the point sources 
contribution, yielding a net increase of the mismatch criterion. 
The 'ideal' case presents an optimistic situation where the ex- 
act contribution of the simulated point sources has been used to 
build the model. This perfect modelling restore the goodness- 
of-fit of the no-ps case. 



point sources increases the variance of the measurement of r. 
In our experiment, the effect is not truly significant (cry shifting 
from 1.84 to 1.91 • 10" 4 ). 

Figure [5] shows the mismatch criterion (from Eq. [16] using 
covariance matrixes binned in €) in the three cases. When no spe- 
cific model of the point source contribution is used, some of their 
emission is nonetheless absorbed by the Smica 'galactic' compo- 
nent, which adjusts itself (via the values of its maximum likeli- 
hood parameters) to represent best the total foreground emis- 
sion. The remaining part is responsible for the increase of the 
mismatch at high I. At the same time, the galactic estimation is 
twisted by the presence of point sources. This slightly increases 
the mismatch on large scales. 

5.2. Galactic foregrounds uncertainties 

We now investigate the impact on the above results of modifying 
somewhat the galactic emission. In particular, we check whether 
a space dependant curvature of the synchrotron spectral index, 
and modifications of the dust angular power spectrum, signifi- 
cantly change the error bars on r obtained in the previous sec- 
tion. 



5.2.1 . Impact of synchrotron curvature 

As mentioned earlier on, the synchrotron emission law may not 
be perfectly described as a single power law per pixel, with 
a constant spectral index across frequencies. Steepening of the 
spectral index is expected in the frequency range of interest. As 
this variation is related to the aging of cosmic rays, it should 
vary on the sky. Hence, the next level of sophistication in mod- 
elling synchrotron emission makes use of a (random) template 
map C(£) to model the curvature of the synchrotron spectral in- 
dex. We then produce simulated synchrotron maps as: 



(3 



ft(f)+ffC(f)log(v/v,) 



(21) 



Betoule et al.: T/S measurements in presence of foregrounds 



13 



Table 4. EPIC-CS measurement for three point sources cases. For the reference case 'no-ps', point sources have neither been 
including in the simulation, nor taken into account in the modelling. The 'ideal' case presents an optimistic situation where the 
exact contribution of the point sources put in the simulation has been used to build the model. The 'no-model' case is a pessimistic 
situation where no effort has been made to model the point sources contribution. 
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^.no-ps j.ideal ^.no-model 


no-py 
07 


^-ideal 




model 


0.001 


1.07 • 10~ 3 1.04 • 10 3 2.00 • 10~ 3 


1.84- 10- 4 


1.91 - 10- 4 


2.49 


• io~ 4 



where a is a free parameter which allows to modulate the ampli- 
tude of the effect (as compared to equation [3}. The right panel 
of figure [6] illustrates the impact of the steepening on the syn- 
chrotron frequency scaling. 

We now investigate whether such a modified synchrotron 
changes the accuracy with which r can be measured. We decide, 
for illustrative purposes, to perform the comparison for EPIC- 
2m, and for r = 0.001 . Everything else, regarding the other emis- 
sions and the foreground model in Smica, remains unchanged. 
Table|5]shows the results of this study in terms of goodness of fit 
and influence on the r estimate. We observe no significant effect, 
which indicates that the foreground emission model of Eq.dTTb is 
flexible enough to accommodate the variation of the synchrotron 
modelling. Even if we cannot test all possible deviation from the 
baseline PSM model, robustness against running of the spectral 
index remains a good indication that results are not overly model 
dependent. 

Table 5. Influence of the running of the synchrotron spectral in- 
dex on component separation in term of goodness of fit and r 
estimates. The study is conducted for the EPIC-2m design, for 
varying amplitude a of the running of the spectral index. No sig- 
nificant variation of the r estimate nor of the likelihood of the 
model is noticed for a remaining in the range allowed by obser- 
vations. 



r cr r 
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r 


-ln.C 







9.78 • 10~ 4 


11.6 


0.001 1.8 • IO -4 


1 


9.62 • 10 4 


11.5 




3 


1.06 • 10~ 3 


11.7 



5.2.2. Level and power spectrum of dust emission 

Similarly, we now vary the model of dust emission and check 
how the main results of section [4] are modified. Measurements 
give some constraints on dust emission on large scales, but 
smaller scales remain mostly unconstrained. Hence, we consider 
here a pessimistic extreme in which we multiply the large scale 
level of the dust by a factor of two, and flatten the power spec- 
trum from a nominal index of -2.5 to -1.9. The power spectra 
corresponding to these two cases are shown in Fig. |6](left panel). 

Running the same component separation pipeline for the 
ground based and the EPIC-2m experiments at their detection 
limit, we find only marginal changes in the measured values of r 
(see table[6}. This result can be interpreted in the following way: 
as the noise of the experiment remains unchanged, the increased 
signal-to-noise ratio allows for a better constraint of the dust pa- 
rameters. Component separation effectiveness depends mainly 



Table 6. Influence of dust polarisation level on component sep- 
aration. We comparison results for a pessimistic (17% intrinsic 
polarisation fraction, flat spectrum) and standard (12% polarisa- 
tion fraction) model of the dust emission. 
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on the coherence of the component, rather than on its overall 
level. 

5.3. Error bar accuracy 

Estimates of the error derived from the FIM (Eq. l20b are ex- 
pected to be meaningful only if the model leading to the likeli- 
hood (Eq. [T4l holds. In particular we assume that processes can 
be modelled as Gaussian. 

We first note that the FIM errors are reasonably compati- 
ble with the difference between input and measured r values, 
which gives confidence that these error estimates are not obvi- 
ously wrong. Nonetheless, we investigate this issue further, us- 
ing Monte-Carlo studies to obtain comparative estimates of er- 
rors, with the EPIC-CS design. Table [7] gives, for two values of 
r and for 100 runs of the Smica pipeline in each case, the aver- 
age recovered value of r, the average error as estimated from the 
Fisher matrix (crJ; ISHER ), and the standard deviation cr^ c of the 
measured values of r. 

For each of the Monte-Carlo runs, a new realization of CMB 
and noise is generated. Simulated galactic foregrounds, however, 
remain unchanged. 

Results show that the FIM approximation give estimates of 
the error in very good agreement with the MC result. Hence, the 
FIM estimate looks good enough for the purpose of the present 
paper, given the number of other potential sources of error and 
the computational cost of Monte-Carlo studies. 

The Monte-Carlo study also allows to investigate the exis- 
tence of a bias. For an input tensor to scalar ratio of 0.01, we 
observe that the measured value of r seems to be systematically 
low, with an average of 9.91 ■ 10 -3 . This we interpret as resulting 
from a slight over-fitting of the data. Still this small bias doesn't 
dominate the error and we are more interested in noise domi- 
nated regime. The overall conclusion of this investigation of er- 
ror bars is that the errors estimated by the FIM are reasonably 
representative of the measurement error. 

5.4. Other cosmological parameters 

The main conclusions of this study are mostly independent of the 
value of all cosmological parameters except t. Within present 
uncertainties indeed, only the value of the reionisation optical 
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Fig. 6. Variations of the galactic foregrounds model. The left panel shows the difference between the default power spectrum of dust 
polarisation B-modes at 150 GHz as modelled by the PSM (solid curve) and a model assuming pessimistic values for the overall 
level and power spectrum index (dotted curve). The right panel shows the dispersion of the synchrotron spectral index for the PSM 
model (in black) and the curved model (in gray). Solid lines present the frequency scaling for the mean values of the spectral index 
and dotted lines for its extremal values. 



Table 7. Monte-Carlo analysis for error bars of the EPIC-CS 
experiment for 2 representative values of r. Sample mean and 
variance are obtained on 100 realizations of noise and CMB. (r) 
denotes the average recovered value of r, (crf ISHER ) the average 
error as estimated from the Fisher matrix, and cr^ the standard 
deviation of the measured values of r 



r 


(r) 






0.01 


9.91 ■ 10- J 
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depth r, which drives the amplitude and p osition of the reionisa - 
tion bump, is critical for our estimations (Colo mbo et al.l l2.Q08). 
Lower r means less accurate measurement of r, and higher r 
better measurement of r. Here we choose a rather conserva- 
tive value of t = 0.07 in agreement with the last measurements 
from WMAP dDunklev et al.ll2008allbh . The value of t, however, 
should affect mainly low resolution and noisy experiments, for 
which most of the information comes from the lowest frequency 
'reionisation' bump in the B-mode spectrum. 

Another issue is that we assume the value of t and n, (and, to 
a less extent, the value of all other cosmological parameters) to 
be perfectly known (setting the shape of the B-mode power spec- 
trum). In fact, uncertainties on all cosmological parameters im- 
ply that the shape will be known only approximately, and within 
a certain framework. Such uncertainties will have to be taken 
into account in the analysis of a real-life data set. Our Smica 
pipeline can be adapted to do this, provided we know the un- 
certainties on the cosmological parameter set. A Monte-Carlo 
approach, in which we assume, for each Smica run, a B-mode 
power spectrum from one of the possible cosmological parame- 
ter sets, will permit to propagate the uncertainties onto the mea- 
surement of r. We expect, however, that this additional error will 
be significantly smaller than that due to the experimental noise. 



6. Conclusion 

In this paper, we presented an investigation of the impact of 
foregrounds on the measurement of the tensor to scalar ratio 
of primordial perturbations. The measurement of r is based on 
the (simulated) observation of the B-mode polarisation of the 
Cosmic Microwave Background by various instruments, either 
in preparation or planned for the future: the Planck space mis- 
sion, a ground-based experiment of the type of Clover, and sev- 
eral versions of a possible dedicated space mission. 

Foreground contamination is modelled and simulated us- 
ing the present development version (v 1.6.4) of the Planck Sky 
Model (PSM). Our main analysis considers the contribution 
from diffuse polarised emission (from the galactic interstellar 
medium modelled as a mixture of synchrotron emission and ther- 
mal emission from dust) and from instrumental noise. The im- 
pact of more complicated galactic foreground emission, and of 
point sources and lensing, is investigated in a second step. 

Our approach uses the Smica component separation method 
on maps of B-modes alone. The method is robust with respect to 
specifics of foreground emission, because it does not rely on an 
accurate representation of foreground properties. That last point 
is demonstrated by varying the input foreground sky, and com- 
paring results obtained with different inputs, without changing 
the analysis pipeline. 

It is shown that for r at the level of r - 0.1, Planck could 
make a meaningful (3(f) detection from B-modes alone. The fi- 
nal sensitivity of Planck for measuring r may be better than what 
is achieved here, as a significant part of the constraining power 
on r should also come from EE/TE for high r. This has not been 
investigated in the present paper, which is more focussed on the 
measurement of low values of r (not achievable with Planck). 
With the various EPIC mission designs, one could achieve de- 
tections at levels of 4-8<x for r = 10" 3 . 

For full-sky, multi-frequency space missions, dealing with 
foregrounds in harmonic space results in a loss of sensitivity by 
a factor 3 to 4, as compared to what would be achievable without 
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foregrounds, even if the covariance of foreground contaminants 
is known. The Smica pipeline allows to achieve performances 
almost as good (within a factor 1 .5), which demonstrates the ef- 
fectiveness of the blind approach, but is still significantly worse 
(factor 3-5) than if there were no foregrounds at all. The loss 
of sensitivity is probably due in part to insufficient localisation 
in pixel space, which results in sub-optimality of the estimator. 
This could (at least in principle) be improved with a localised 
processing. 

For the most ambitious EPIC space mission, we find that 
our main conclusions are not modified significantly when taking 
into account the contamination of primordial B-modes by extra- 
galactic point sources, by gravitational lensing, or when simu- 
lating a more complicated galactic emission. In contrast, we find 
that the measurement of r from the ground with few frequency 
channels can be severely compromised by foregrounds, even in 
clean sky regions. 

The joint analysis of such ground-based data together with 
those from less sensitive experiments covering a wider fre- 
quency range, as the Planck data, permits to improve the con- 
straints on r. Still, the result from a combined analysis of Planck 
and of a small patch observed from the ground at few frequencies 
cannot match what is obtained using sensitive measurements on 
the whole frequency range. 

This makes a strong case for sensitive multi-frequency ob- 
servations, and thus probably also for a space mission, as obser- 
vations from the ground are severely limited (in frequency cov- 
erage) by atmospheric absorption and emission. This conclusion 
is further supported by the fact that a space mission mapping the 
same clean region (about 1 % of the sky), but with the full fre- 
quency range allowed by the absence of atmosphere, makes it 
possible to deal with diffuse foregrounds very efficiently. 

Such a deep field mission would, in that respect, outperform 
a comparable full-sky experiment. The results obtained in the 
present study, however, do not permit to conclude whether a full 
sky or a deep field mission would ultimately perform better. A 
strategy based on the observation of a small patch seems to offer 
better prospects for measuring r with an harmonic-space based 
version of Smica, but also seems to be more impacted by small 
scale contamination than all-sky experiments, and is in particular 
quite sensitive to the lensing effect. Further developments of the 
component separation pipeline could improve the processing of 
both types of datasets. 

As a final comment, we would like to emphasise that the 
present study is, to our knowledge, the first one which de- 
signs, implements effectively, and tests thoroughly on numer- 
ous simulations a component separation method for measur- 
ing r with CMB B-modes without relying much on a physical 
model of foreground emission. The method is shown to be ro- 
bust against complicated foregrounds (pixel-dependent and run- 
ning synchrotron spectral index, multi-template dust emission, 
polarised point sources and lensing). It is also shown to pro- 
vide reliable errors bars on r by comparing analytical error bars 
(from the FIM) to estimates obtained from Monte-Carlo simula- 
tions. Although more work is needed for the optimal design of 
the next B-mode experiment, our results demonstrate that fore- 
grounds can be handled quite effectively, making possible the 
measurement of r down to values of 0.001 or better, at the 5-6cr 
level. 

Certainly, next steps will require fully taking into account 
small scale contaminants, partial sky coverage effects, and prob- 
ably some instrumental effects in addition to diffuse foregrounds. 
For this level of detail, however, it would be mandatory to refine 
as well the diffuse foreground model, using upcoming sensitive 



observations of the sky in the frequency range of interest and 
on both large and small scales. Such data will become available 
soon with the forthcoming Planck mission. 
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Appendix A: Parameterisation the foreground 
component and choice of a mask 

In this appendix, we discuss in more detail the dimension D of 
matrix used to represent the covariance of the total galactic emis- 
sion, and the choice of a mask to hide regions of strong galactic 
emission for the estimation of r with Smica. 

A. 1 . Dimension D of the foreground component 

First, we explain on a few examples the mechanisms which set 
the rank of the foreground covariance matrix, to give an intu- 
itive understanding of how the dimension D of the foregrounds 
component used in Smica to obtain a good model of the data. 

Let's consider the case of a 'perfectly coherent' physical pro- 
cess, for which the total emission, as a function of sky direction 
£ and frequency v, is well described by a spatial template multi- 
plied by a pixel-independent power law frequency scaling: 

S Y (f) = So(f){£j (A.l) 

The covariance matrix of this foreground will be of rank one 

and R s = [AA t var(So)], with Aj = (jf-f- Now, if the spectral 
index B fluctuates on the sky, B{%) - B + 6B(^), to first order, the 
emission at frequency v around vo can be written: 

S v (0 * Sotf) f-f + S (& f-f TO f (A.2) 
\vo/ \vo/ \ v / 

This is not necessarily the best linear approximation of the 
emission, but supposing it holds, the covariance matrix of the 
foreground will be of rank two (as the sum of two correlated 
rank 1 processes). If the noise level is sufficiently low, the varia- 
tion introduced by the first order term of Eq. IA.2l becomes truly 
significant, we can't model the emission by a mono-dimensional 
component as in Eq lA.ll 

In this work, we consider two processes, synchrotron and 
dust, which are expected to be correlated (at least by the galactic 
magnetic field and the general shape of the galaxy). Moreover, 
significant spatial variation of their emission law arises (due to 
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Table A.l. Bayesian information criterion of 3 models with in- 
creasing dimension of the galactic component for the EPIC-2m 
mission. The selected value D — 4 correspond to a minimum of 
this criterion. 



cosmic aging, dust temperature variation ...), which makes their 
emission only partially coherent from one channel to another. 
Consequently, we expect that the required dimension D of the 
galactic foreground component will be at least 4 as soon as the 
noise level of the instrument is low enough. 

The selection of the model can also be made on the basis of a 
statistical criterion. For example, Table |A.f1 shows the Bayesian 
information criterion (BIC) in the case of the EPIC-2m exper- 
iment (r = 0.01) for 3 consecutive values of D. The BIC is a 
decreasing function of the likelihood and of the number of pa- 
rameter. Hence, lower BIC implies either fewer explanatory vari- 
ables, better fit, or both. In our case the criterion reads: 

BIC = -21nX + Jtln^ w q 

i 

where k is the number of estimated parameters and w q the effec- 
tive number of modes in bin q. Taking into account the redun- 
dancy in the parameterisation, the actual number of free parame- 
ters in the model is 1 + F x D + QD(D + l)/2 - D 2 . However, we 
usually prefer to rely on the inspection of the mismatch in every 
bin of {, as some frequency specific features may be diluted in 
the global mismatch. 

A.2. Masking influence 

The noise level and the scanning strategy remaining fixed in the 
full-sky experiments, a larger coverage gives more information 
and should result in tighter constraints on both foreground and 
CMB. In practice, it is only the case up to a certain point, due to 
the non stationarity of the foreground emission. In the galactic 
plane, the emission is too strong and too complex to fit in the 
proposed model, and this region must be discarded to avoid con- 
tamination of the results. The main points governing the choice 
of an appropriate mask are the following: 

- The covariance of the total galactic emission (synchrotron 
and dust polarised emissions), because of the variation of 
emission laws as a function of the direction on the sky, is 
never exactly modelled by a rank D matrix. However it is sat- 
isfactorily modelled in this way if the difference between the 
actual second order statistics of the foregrounds, and those 
of the rank D matrix model, are indistinguishable because of 
the noise level (or because of cosmic variance in the empiri- 
cal statistics). The deviation from the model is more obvious 
in regions of strong galactic emission, hence the need for a 
galactic mask. The higher the noise, the smaller the required 
mask. 

- Smica provides a built-in measure of the adequacy of the 
model, which is the value of the spectral mismatch. If too 
high, the model under-fits the data, and the dimension of the 
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Table A.2. Estimation of the tensor to scalar ratio with two dif- 
ferent galactic cuts in the EPIC-2m experiment. 



foreground model (or the size of the mask) should be in- 
creased. If too low, the model over-fits the data, and D should 
be decreased. 

- Near full sky coverage is better for measuring adequately the 
reionisation bump. 

- The dimension of the foreground component must be smaller 
than the number of channels. 

If the error variance is always dominated by noise and cos- 
mic variance, the issue is solved: one should select the smaller 
mask that gives a good fit between the model and the data to 
minimise the mean squared error and keep the estimator unbi- 
ased. 

If, on the other hand, the error seems dominated by the con- 
tribution of foregrounds, which is, for example, the case of the 
EPIC-2m experiment for r - 0.001, the tradeoff is unclear and 
it may happen that a better estimator is obtained with a stronger 
masking of the foreground contamination. We found that it is 
not the case. Table [AT2l illustrates the case of the EPIC-2m ex- 
periment with the galactic cut used in Sect. [4] and a bigger cut. 
Although the reduction of sensitivity is slower in presence of 
foreground than for the noise dominated case, the smaller mask 
still give the better results. 

We may also recall that the expression (0 of the likelihood 
is an approximation for partial sky coverage. The scheme pre- 
sented here thus may not give fully reliable results when mask- 
ing effects become important. 

Appendix B: Spectral mismatch 

Computed for each bin q of I, the mismatch criterion, 
w q K(VK q , R ? (#*)), between the best-fit model R ? (#*) at the point 

of convergence 9*, and the data R fy , gives a picture of the good- 
ness of fit as a function of the scale. Black curves in Figs. IB. 11 
and IB. 21 show the mismatch criterion of the best fits for Planck 
and EPIC designs respectively. When the model holds, the value 
of the mismatch is expected to be around the number of degrees 
of freedom (horizontal black lines in the figures). We can also 
compute the mismatch for a model in which we discard the CMB 
contribution w q K(R q , R q (9*) - R™ B (r*)). Gray curves in Figs. 
IB. 11 and IB. 21 show the mismatch for this modified model. The 
difference between the two curves illustrates the 'weight' of the 
CMB component in the fit, as a function of scale. 

Figure IBTTI shows the results for Planck for r = 0.3 and 0.1. 
The curves of the difference plotted in inclusion illustrate the 
predominance of the reionisation bump. In Fig. IB. 21 we plot the 
difference curve on the bottom panels for the three experiments 
for r = 0.01 and r = 0.001. They illustrate clearly the differ- 
ence of sensitivity to the peak between the EPIC-LC design and 
the higher resolution experiments. In general it can be seen that 
no significant contribution to the CMB is coming from scales 
smaller than I = 150. 
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Fig. B.l. Those plots present the distribution in i of the mis- 
match criterion between the model and the data for two values 
of r for Planck. On the grey curve, the mismatch has been com- 
puted discarding the CMB contribution from the Smica model. 
The difference between the two curves, plotted in inclusion, il- 
lustrates somehow the importance of the CMB contribution to 
the signal. 
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Fig. B.2. Mismatch criterion for r = 0.01 (top) and r = 0.001 (bottom). In each plot, the top panel shows the mismatch criterion 
between the best fit model and the data (black curve) and the best fit model deprived from the CMB contribution and the data (gray 
curve). Solid and dashed horizontal lines show respectively the mismatch expectation and 2 times the mismatch expectation. The 
difference between the gray and the black curve is plotted in the bottom panel and gives an idea of the significance of the CMB 
signal in each bin of I. 



